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EFFECT OF SURFACE ROUGHNESS ON I- YDRODYNAMIC BEARINGS 


Bankim C. Majumdar and Bernard J. Hamrock 
National Aeronautics and Space Aoministration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

A theoretical analysis on the performance of hydrodynamic oil bearings 
is made considering surface roughness effect. The hydrodynamic as well as 
asperity contact load is found. Assuming the surface height distribution as 
Gaussian the contact pressure is calculated. The average Reynolds equation 
of partially lubricated surface is used to calculate hydrodynamic load. An 
analytical expression for average gap is found and is introduced to modify 
the average Reynolds equation. The resulting boundary value problem is then 
solved numerically by finite difference methods using the method of succes- 
sive over-relaxation. The pressure distribution and hydrodynamic load capa- 
city of plane slider and journal bearings are calculated for various design 
data. The effects of attitude and roughness of surface on the bearing per- 
formance are shown. The results are compared with similar available solu- 
tion of rough surface bearings. It is shown that (1) the contribution of 
contact load is not significant and (2) the hydrodynamic and contact load 
increase with surface roughness. 


INTRODUCTION 

The classical theory of hydrodynamic lubrication given by Reynolds does 
not consider the surface roughness of the elements having relative motion. 
This theory is applicable when the bearing surfaces are completely separated 
by a thick lubricant film. It is known that the bearing load supporting 
ability increases with decrease in the lubricant film. When the load is 
very high and the film thickness is small there is a possibility of asperity 
contact. The method of computation of surface contact load is extremely 
complex. A few theories concerning the contact of nominally flat surfaces 
are available [1-4]. In addition to this contact load, the lubricant film 
between two surfaces having relative motion will also carry a load. This 
load can be called the hydrodynamic load of a rough surface bearing. There 
exists two main approaches [5-7] for computation of the hydrodynamic load. 
Patir and Cheng [5,6] used the flow simulation method of a randomly genera- 
ted rough surface with known statistical properties over the bearing area. 
Tender [7] studied the lubrication of a rough surface by a Monte Carlo 
method. 

The present study uses a flow model S'imilar to Patir and Chenq [5]. 

This study differs from [6] so far as calculation of average gap ny is 
concerned. This gap height is obtained analytically in terms of nominal 
film thickness, h. It is then introduced in the basic equation to modify 
average Reynolds equation. Moreover, the computation of contact load in 
addition to hydrodynamic load is included. The partial differential equa- 
tion is solved numerically by finite difference methods satisfying the ap- 
propriate bearing boundary conditions. The effect of rc ghness parameter, 
attitude, and ratio of correlation length of asperity on the hydrodynamic 
load capacity is shown. 
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NOTATIONS 

Ap apparent area of contact [m^] 

B width of plane slider [m] 

C nominal clearance in journal bearing [m] 

^a’^b modulus of elasticity of material of surfaces a and b LN/m^] 



e eccentricity [mj 

h,H nominal film thickness [m], H = h/h^ for slider bearing, 

H = h/C for journal bearing 

hj^,h 2 nominal maximum and minimum film thicknesses of slider bearing [mj 
h^ local film thickness [mJ 

h^ average film thickness [m] 

K a constant 

L length of bearing [m] 

N number of asperities per unit area 

n attitude, n = h^/n^ 

p,P mean hydrodynamic pressure [N/m'^j, P * (ph^j/IbnuB) for 

slider bearing, P = (pC^)/(6r,uk) for journal bearing 
contact pressure [N/m ] 

R journal radius [m] 

t time [s] 

Ua.Up.u velocity of surfaces [m/sj, = u, = 0 
v,V squeeze velocity,v = L-(3h)/(n)] [m/sj, V = (Bvj/lh^u) tor 

slider bearing, V = (Rv)/(Cu) for journal bearing 
^cj^c asperity contact loud [N], = w^/LBE' for slider bearing, 

^c = w^/2LRE' for journal bearing 
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? 7 

^h’^h hydrodynamic load [N], = (w^hp/(6nLB u) for slider bearing, 

2 2 

for journal bearing, * (w^C )/(6oLR u) for journal 
beari ng 

x,y,X,Y coordinates [m], X = x/B, Y = y/L 

e mean radius of curvature of asperities [m] 

Y ratio of x and y correlation length (surface pattern parameter) 

A roughness parameter, a = h^/o for slider bearing, A = C/o for 

journal bearing 

«a,«b»® roughness amplitudes of surfaces measured from their mean levels 
[m], 6 = combined roughness, ® 65 

t eccentricity ratio, e = e/C 

2 

n coefficient of absolute viscosity of lubricant i.Ns/in ] 


9,62 


Poisson's ratio of material of surfaces a and b 


angular coordinate (rad), e = x/R, «2 = angular coordinate where 


f i Im breaks 

standard deviation of combined roughness 6, 0 
composite rms roughness 



0^,05 standard deviations of roughness functions 6^ and 6^ 
4 ' attitude angle [rad] 


^ pressure flow factors 

X y 

shear flow factor 

THEORY 


The average Reynolds equation for partial hydrodynamic lubrication is 
given by Patir and Cheng [b] and it can be written for the surfaces shown in 
figure 1 as 


3 __ 

3X 


h^ ap\ 
X j 



tr_ ^ 

'y 12n 3 y 


U, + U, 3Kt U, - U, 3^1^ 3h-r 

b_T^_a b^__s^__T 

2 3x 2 3x at ' 


where 
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a is the standard deviation 

i>x and 0y are pressure flow factors 

^5 is shear flow factor 

The ave.'-age gap is calculated from 


\ = (h + 6)f(6)d« 
' <-h 


( 2 ) 


where f(«) is probability density function of combined roughness, 6. The 
flow factors 0^, 0y will approach 1 as h/o approaches whereas 05 
will be equal to zero for a large value of h/o. The average gap 
is seen to be functions of combined ghness and the probability density 
function of 6. The study to foil w nil deal with when both surfaces have 
the same roughness structure and same rms roughness. 

To obtain a solution of equation (1) for a particular bearing configur- 
ation one has to find 0^, 0y, 0$, and Tij beforehand. The 0^, 0y, and 
05 are dependent on the roughness geometry of the bearing surface."^ The 
flow factors not only depend on h/o as mentioned above, they are functions 
of the statistical properties, such as the frequency density of roughness 
heights and the directional properties of the asperities. The height dis- 
tribution is assumed to be Gaussian. The flow factors 0x and 0y are 
obtained by Patir and Cheng [5] through flow simulation of a rough surface 
having Gaussian distribution of surface height. These are used in the pres- 
ent case. However, the average gap hj is calculated in the following 


way: 

For a Gaussian distribution the normal probability density function of 

6 i s 


f(6) 




(3) 


where o is tl.e standard deviation. 

Substituting equation (3) into equation (2), 



After performing the above integration, we get 


hj as 


TT h 
hf = 2 


1 + erf 



0 



e 


(4) 
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SI nee 


/ *> 2 

e”'" dr • i[l - erf(h)] 


Differentiating hj with respect to x, one obtains 

aTi, 


'T 1 

3X ■ 7 


1 + erf 


fe) 


ill 

3X 


Similarly 


aTi-^ 

TT 


- f -"fe) 


3h 

3t 


The flow simulation factors Px» ^y» ^s given by [5,6J, 

for Y < 1 


*x = 1 - Cje 


-rl- 
0 


(5) 


( 6 ) 


and 


*x = ‘ * Cidf for , > 1 


,, “1 a,n‘ 


for - < 5 
0 — 


h = '' 2 ® 


-“ 4 O 


for - > 5 

Cl 


^(^ '') ' ^(^ 7 ) 


(7) 


where Ci, r, Ai, A 2 , ai, 02 , 03 , and 04 are constants and tabulated in 
[ 6 ], and y Is defined as ratio of lengths at which the auto-correlation 
function of the x and y profiles reduce to 50 percent of the initial 
value. This y can oe thought of as the length-to-width ratio of a repre- 
sentative asperity. For isotropic roughness y = 1* 

Having known the flow factors and Tii-; and assuming u^ = 0 and 
Ug = u equation ( 1 ) for constant n and when both surfaces have same 
roughness structure can be written as 
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It has been found in [faj that for the type of model as mentioned earlier 
{i.e., when both surfaces have same roughness structure) there is no 
effect. 

Let us now attempt to find the solution of equation (8) for the infi- 
nitely long plane slider, finite plane slider, and infinitely long journal 
beari ng. 


Infinitely Long Plane Slider Bearing 


Figure Z shows a slider bearing. If length of the slider is very long 
compared to B, 8o/3y = 0. Hence equation (8) can be written as 


3X 




Using the following substitutions: 

h 


ph: 


y ^ . 2 u h ^ , 0 *'"2 

X = g. A = g , H = 1 ^^, and P = ^ 


(9) 


equation (9) becomes 


It) ‘ i 


1 + erf 



aX “ ho 




(10) 


SI nee V = -( 3h/ at) . 

Fxpanding equation (10), 


* 


H aX* 7 


• "'(^ij|s - ■> 


Here 


V = 


Bv 

h^u 


( 11 ) 


The dimensionless parameter V can be visualized as squeeze effect and it 
comes from t.he flow continuity. 

For a slider bearing, the dimensionless film thickness, H is given by 

H = n - (n - 1 )X (Iz) 

rtPei e n = h]^/h 2 . 

The bearing boundary conditions are 
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and 


P = 0 at X = 0 
P * 0 at X = 1 


(13) 


Finite Slider Bearing 

The dimensionless differential equation in this case will be 



where Y = y/L. 

The film thickness H, and the bearing boundary conditions are still 
given by equations (12) and (13) with additional boundary conditions at the 
sides. 


I nf initely Long Journal Bearing 

A journal bearing as shown in figure 3 when rotating with velocity u 
carries a load w. If the length of the bearing is large compared to other 
dimensions, then there will be no side leakage. Thus when L/R is large, 
3p/3y = 0. Therefore the basic equation will be equation {'i). This can be 
nondimensional ized with the following substitutions: 

2 

X U h . C , n pC 
G-p,H-^, A=^, and P - 


The dimensionless equation is 



(lb) 


where 

H = 1 + E cos ©, £ = ^ ~ 

The boundary conditions are; 

P = 0 at e = 0 
P . |£ . 0 at 0=0^ 


( 16 ) 


( 17 ) 
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where is the angular coordinate at which the film breaks. This bound- 
ary condition is known as Reynolds condition. Using the expression of H, 
equation (15) can be reduced to 


2 

<)^(1 + € cos e)^ ^ - 3c0^(l + C cos ®)^ sin e 
96 


3 aP 

(1 ♦ c cos 


1 


= ? 





sin e - V) 


(lb) 


Equations (11), (14), and (18) are solved numerically by finite differ- 
ence methods with successive over-relaxation factor satisfying the appropri- 
ate boundary conditions for infinitely long plane slider, finite plane 
slider, and infinitely long journal bearing, respectively. 

Calculation of Hydrodynamic Load 

With the known hydrodynamic pressure distribution the hydrodynamic load 
can be calculated from: 

For infinitely long plane slider - 

“h ' L p dx 


or 


where 


faoLB^u 



P dX 




^ 

bnLB U 


For finite plane slider - 


"h 



(ly) 


(^ 0 ) 


v.here w^ is defined the same as for the infinite slider. 

For infinitely long journal bearing the two components of loads are - 



y 


w 


hr 



p cos e cle 


and 


'h« 



p sin e de 


These can be written in dimensionless form as 


w 


hr 



P cos e de 


and 




sin e do 


( 21 ) 


The total hydrodynamic load is 


where 





" ‘ 1 1.2 
faoLR u 


The attitude angle 



( 22 ) 


(23) 


The integrations of equations (19), (20), and (21) are performed numer- 
ically by Simpson's rule. 

Calculation of Asperity Contact Load 

Using a Gaussian distribution of csperity height, tne contact ioaa can 
be evaluated from the nominal contact pressure. When both surfaces are 
rough, the nominal pressure can be obtained from Greenwood and Tripp [4j. 
This is written as 


^c ■ ^ ^ 


(24) 


where 



lu 



N number of asperities per unit area 

a mean radius of curvature of asperities 

E' composite modulus of elasticity 

and 


2 

C 1 i -{A 

Ft/p = / (s - A ) ' e ' ' 'ds 

A 


The value of K varies from 0.003 to 0.0003 for the range of o/a 
between 0.01 and 0.001. The function Fij/2 has been calculated in [4] 
and is reproduced below for various A. 

A 

0 

.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

The contact load is 


^ 5/2 

0.61664 
. 24040 
.08056 
.02286 
.00542 
.00106 
.00017 
.00002 
.00000 


where Ap is the apparent area of contact. 
Therefore for plane slider bearing 


or 


w, = Lbp. 
c c 


«c = 


^ = lET 


[cb 


where 
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similarly for journal bearing W^. is given by 


70 ^ 


where W^. ■ 2LRpc« 

The above two dimensionless loads and are defined with 
different physical parameters. Therefore these cannot be simply added in 
this form. In the following section an example is taken to show tiie effect 
of and Wh« The discussion to follow is concerning W^. 

RESULTS AND DISCUSSION 

The hydrodynamic load has been calculated using the foregoing method 
for infinitely long plane slider, finite slider, and infinitely long journal 
bearings. The example given later shows that the asperity contact load is 
much smaller than the hydrodynamic load. Therefore the results are given 
with respect to only. 

Figures 4 to 10 show the results of plane slider bearings and the re- 
sults for the journal bearing are given in table I and figure 11. It ma> be 
noted the results corresponding to A * 6 approach the smooth bearing so- 
lutions. 


Plane SI ider Bearing 

The effect of various parameters, namely, n, V, L/B, A and > on 
Wh is d iscussed in the following paragraphs. 

(1) Effect of roughness parameter a: Each of the figures shows that 

the hydrodynamic load is increased with decrease in the A. This increase 
is due to pressure flow effect and py) and the extra term in the 
average gap The hydrodynamic load also increases very sharply at 

small A. 

(2) Effect of attitude n: when local squeeze term V is neglected, 

the load capacity is maximum at n = 2 for most values of A . The varia- 
tion of load with n for a particular a is more or less similar to that 
of a smooth bearing. It may be mentioned that an infinitely long plane 
slider having smooth surfaces carries a maximum load when n . 2,2. In the 
present case the load reaches a maximum value when n varies between 2.0 
and 2.b. 

(3) Effect of local squeeze velocity V: The effect of squeeze veloc- 

ity is always to enhance tne load. In rough surface bearings (tigs. 4 

and b) the similar effect is also seen. From figure o it may be noticed 
thai a bearing operating on small attitude gives high load. This indicates 
tnat for a bearing having nearly parallel film, the squeeze velocity plays a 
significant role so far as load capacity is concerned, 

(4) Effect of length to width ratio L/B: The hydrodynamic load ifi- 

creases with the L/B ratio (fig. y). This is expected from the physical 
point of view. However, load does not : 'is i stent ly change with the cnange 
of L/B ^ .'tio. 

(5) Ei'fect of surface pattern parameter The i parameter repre- 
sents tne directional properties of surfaces. As y is tne ratio of 
length-to-width of an asperity and as the orientation of asperity disturbs 
the flow, its effect on the load .s seen to predominant (fig. lO). As per 
definition purely transverse, isotropic and longitudinal rougnness patterns 
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correspond to y » 0, 1 and «, respectively. As the values of pressure 
flow factor are known from 1/6 to 6, these are used here. These two 
values adequately represent transverse and longitudinal patterns. For a 
square bearing the maximum and minimum load are obtained with isotropic and 
longitudinal surface roughness, respectively. Whereas the surfaces with 
transverse pattern give intermediate values. 

Approximate Equation 

From a near optimum slider (n s 2) it has been found that the load ca- 
pacity of isotropic rough and smooth surfaces can be closely approximated by 
the following simple relationship: 


h( smooth) 
h(rough) “ ^ _ o.77 


The maximum error introduced by this approximation is about 5 percent. 

Journal Bearing 

The hydrodynamic load, attitude angle, and the location of lubricating 
film breakdown for an infinitely long journal bearing using Reynolds bour.d- 
ary conditions are given in table I and figure 11. The load capacity in- 
creases with decrease in the C/o (or A ) ratio, but the extents of lubrica- 
ting film in most cases do not change. The attitude angle drops very slowly 
with increase in the roughness. Due to this reason the variation of 4 > 
and ©2 given in table I. The variation of load, attitude angle, and 
the extent of film with eccentricity ratio is very much similar to that of a 
smooth bearing. 


Calculation of Total Load (w^ + w^-) 

The dimensionless hydrodynamic and contact loads are defined in such a 
way one cannot adc them directly. Ihe relative load sharing ability is 
shown with an illustrative example. 

E X atiif) I e : 

A plane slider bearing made of steel is operating under the following 
cond 1 t ions: 

Length of bearing L = bU mm 

width of bearing B = bO mm 

Minimum film thickness h^' = bxlO“*^ m 

Att itude n = 2.U 

Sliding speed u = b m/s 

Absolute viscosity of oil n = 0.20 Ns/m^ 

The hydrodynamic and contact load for the above bearings are calculated. 
Assume [' = 2.2xlol^ N/';r/, 

The hydrodynamic load is given by the expression: 
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w. 


bnLB^u 


Wj,LNj 


Substituting the above data in this expression 

= 30x10^ W^[N] 

The contact load can be written as 

= LBKE'F^,2tN] 

Although will be slightly different from h^ to h]^, F5/2 

is used here. Taking K = 0.003, assuming isotropic surface roughness and 

using the above data 


w^= 165x10^ 

For various roughness parameters the hydrodynamic and contact load are 
calculated for isotropic surface roughness and are shown below: 


Al^h^/o) 

w^[kN] 

w^LkN] 

w /w, 
c h 

1 

453 

133 

0.294 

2 

360 

8.94 

.0248 

3 

321 

.28 

.00081 

6 

300 

0 

0 


comparison of Results 

The hydrodynamic load capacity for the square plane slider with n = 2 
and Y = 1 is compareo with that of Patir and Cheng [6] result. These are 
shown in table II. The Patir and Cheng results are slightly higher than 
those obtained from the present method of solution. All the data in 
table II are in agreement within 20 percent. The load capacity correspond- 
ing to other values of A are consistently higher. This may be due to 
assumption of poj_ynomial density function for the Gaussian function in the 
calculation of hj in [6]. T_he present calculations are, however, made 
from tne exact expression of hj (eq. (4)). Again the numerical calcula- 
tion are done with an accuracy of 0.01 percent of the difference of integra- 
ted pressures of two successive iterations. Therefore, it is believed that 
the present data give good accuracy, 

CONCLUSIONS 

Loads due to asperity contact as well as hydrodynamic pressure were 
considered in studying the effect of surface roughness on hydrodynamic bear- 
ings. The height distribution was assumed to be Gaussian aiid the coritact 
pressure was calculated. The average Reynolds equation with flow simulation 


u 


of a randomly generated rough surface was used to calculate the hydrodynamic 
load. An analytical expression for the average gap was found and was intro- 
duced in the modified average Reynolds equation. The pressure distribution 
and hydrodynamic load capacity of plane slider and journal bearii^s are es- 
timated for various design data. The effect of attitude and ratio of film 
thickness to standard deviation of a surface was shown. The following con- 
clusions are drawn from the above analysis. 

1. Both the hydrodynamic and contact load increase with increase in the 
surface roughness. 

2. For a square plane slider maximum hydrodynamic load is observed with 
isotropic surface roughness. 

3. The hydrodynamic load approaches to a smooth bearing solution 
when A = 6 and the contact load approaches zero when A = 4. 

4. Although the hydrodynamic load of journal bearing varies with sur- 
face roughness, there is little variation of attitude angle and the location 
of the point of film breakdown. 

5. The approximate equation proposed for plane slider bearing having 
n = 2 may be used for design calculations. 
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TABLE II. - COMPARISON OF HYDRODYNAMIC LOAD CAPACITY OF 


FINITE SLIDER BEARING (L/B = 1, n = 2, y = 1) 


m 

Wh 





0.0132 


2 

.0146 


3 

.0130 


6 

.0115 










HYDRODYNAMIC LOAD. 


Figure 3. - A journal bearing. 



SURFACE ROUGHNESS PARAMETER, A 

Figure 4. - Variation of hydrodynamic load with surface roughness 
parameter for various attitudes. 




SURFACE ROUGHNESS 


PARAMETER. A 

Figure 5. - Variation of hydrodynamic 
load with surface roughness parameter 
for various attitudes. 


HYDRODYNAMIC LOAD. 
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SURFACE ROUGHNESS PARAMETER. A 

Figure 6. - Variation of hydrodynamic load with surface roughness 
parameter for various attitudes. 



HYDRODYNAMIC LOAD. 



SURFACE ROUGHNESS PARAMETER, A 

Figure 7. - Variation of hydrodynamic load with surface roughness 
parameter for various attitudes. 


HYDRODYNAMIC LOAD. 



Figure 8. - Variation of hydrodynamic load with surface rough- 
ness parameter for various attitudes. 




HYDRODYNAMIC LOAD. 



SURFACE ROUGHNESS PARAME7IR, A 

Figure 10. - Variation of hydrodynamic load with rough- 
ness parameter for various surface patterns. 



Figure 11. - Variation of load of an infinitely long journal bearing 
with roughness parameter for various eccentricity ratios (V » 0). 
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